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" B schreibung 

[0001] Verfahren zum modellgestutzten Bestimmen der in die Zy Under einer Brennkraftmaschine einstromenden 
Frischluftmasse bei externer Abgasruckfuhrung 

[0002] Die Erfindung betrifft ein Verfahren zum modellgestutzten Bestimmen der in die Zylinder einer Brennkraftma- 
schine einstromenden Frischluftmasse bei externer Abgasruckfuhrung gemaS dem Oberbegrifl von Patentanspruch 1 . 
[0003] Motorsteuerungssysteme fur Brennkraftmaschinen, die mit Kraftstoffeinspritzung arbeiten, verwenden als 
MaB fur die Motorlast entweder die vom Motor angesaugte Luftmasse oder den Saugrohrdruck. 
[0004] Luftmassengef uhrte Motorsteuerungssysteme benotigen also fur die Berechnung der Einspritzzeit diejenige 
Frischluftmasse, die aus dem Saugrohr heraus in die jeweiligen Zylinder stromt. 1st die Brennkraftmaschine mit einer 
externen Abgasruckf uhrvorrichtung ausgestattet, mit deren Hilfeein Teildes Abgases in den Ansaugtraktzuruckgefuhrt 
wird s so stelit im Instationarbetrieb, der z.B. durch eine Anderung der externen Abgasruckfuhrrate und/oder der Dros- 
selklappenstellunghervorgerufen werden kann, das Signal eines stromauf warts der Drosselklappe angeordneten Luft- 
massenmessers kein MaB fur die tatsachliche Fullung der Zylinder dar. Die Fullung andert sich entsprechend der 
Anderung des Frischgaspartialdruckes im Saugrohr, so daB zur exakten Bestimmung der LastgroBe der Absolutdruck 
und die RartialdrOcke im Saugrohr bekannt sein mussen. 

[0005] Saugrohrdruckgef uhrte Motorsteuerungssysteme konnen uber eine Druckmessung nur den Absolutdruck im 
Saugrohr erfassen. Da die Kraftstoffzumessung jedoch durch den Frischgaspartialdruck bestimmt wird, mussen der 
Frischgas- und bei Abgasruckf uhrung auch der Restgaspartialdruck auf geeignete Weise bebestimmt werden. 
[0006] Auf physikalischen Ansatzen basierende modellgestutzte. Berechnungsmethoden stellen einen guten Aus- 
gangspunkt zur genauen Bestimmung der Luftmasse dar. 

[0007] In der nicht vorveroffentlichten deutschen Patentanmeldung P 195 13 601.2 derselben Anmelderin ist ein 
Verfahren zum modellgestutzten Bestimmen der in die Zylinder einer Brennkraftmaschine ohne Abgasruckf Ohrung 
einstromenden Luftmasse beschrieben. Mit Hilfe eines Saugrohrfullungsmodells werden dabei die GroBen Saugrohr- 
druck, Luftmassenstrom an der Drosselklappe und Luftmassenstrom in den Zylinder aus dem Drosselklappenoffnungs- 
grad und der Motordrehzahl angegeben und eine LastgroBe berechnet, auf dessen Grundlage die Einspritzzeit be- 
stimmt wird. 

[0008] Auf der Basis des Signals des Lasterfassungssensors, d.h. des Luftmassenmessers bzw. des Saugrohrdruck- 
sensors, werden uber einen Modellregelkreis die Wirkung von Parameterfehlern des Saugrohrfullungsmodells redu- 
ziert, indem die Abweichungen zwischen der MeBgroBe und der entsprechenden ModellgroBe minimiert werden. Die- 
ser Zugang gestattet eine permanente Modellkorrektur im stationaren und instationaren Motorbetrieb, so daB die Last- 
groBe der Brennkraftmaschine mit hoher Genauigkeit geschatzt werden kann. 

[0009] Aus der DE 39 19 488 C2 ist eine Vorrichtung zur Regelung und zur Vorausbestimmung der Ansaugluftmenge 
einer saugrohrdruckgefuhrten Brennkraftmaschine bekannt, bei der der Drosselklappenoffnungsgrad und die Motor- 
drehzahl als Grundlage zur Berechnung des derzeitigen Wertes der in den Brennraum der Maschine eingesaugten 
Luft verwendet werden. Diese berechnete, gegenwartige Ansaugluftmenge wird dann als Grundlage zur Berechnung 
des vorausbestimmten Wertes fur die Ansaugluftmenge, die in den Brennraum der Maschine zu einer bestimmten Zeit 
von dem Punkt an, an dem die Berechnung ausgefuhrt wurde, einzusaugen ist, benutzt. Das Drucksignal, das strom- 
abwarts der Drosselklappe gemessen wird, wird mit Hilfe von theoretischen Beziehungen korrigiert, so daB eine Ver- 
besserung der Bestimmung der angesaugten Luftmasse erreicht und damit eine genauere Berechnung der Einspritz- 
zeit mogiich ist. 

[0010] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren anzugeben mit dem bei einer Brennkraftmaschine 
mit Abgasruckfuhrung der Frisch- und Restgaspartialdruck im Saugrohr im stationaren und instationaren Motorbetrieb 
moglichst exakt bestimmt und als Grundlage fur eine genaue Berechnung der Einspritzzeit herangezogen werden kann. 
[0011] Diese Aufgabe wird gemaG den Merkmalen des Patentanspruches 1 gelost. Vorteilhafte Weiterbildungen 
finden sich in den Unteranspruchen. 

[001 2] Ein Ausf uhrungsbeispiel des erf indungsgemaBen Verfahrens wird anhand der nachfolgenden schematischen 
Zeichnungen beschrieben. Dabei zeigen: 

Figur 1 eine Prinzipskizze zum Saugsystem einer Otto-Brennkraftmaschine einschlieBlich deT entsprechenden 

Modell- und MeBgroBen, 
Figur 2 die DurchfluBfunktion und die dazugehorige Polygonzugapproximation, 

Figur 3 eine Prinzipdarstellung eines Modellregelkreises zur Bestimmung der in die Zylinder der Brennkraftma- 
schine einstromenden Frischluftmasse bei aktiver Abgasruckfuhrung und 

Figur 4 eine Prinzipdarstellung eines Modellregelkreis zur Korrektur des Zusammenhanges zwischen Saugrohr- 
druck und Luftmasse wahrend abgeschalteter Abgasruckfuhrung 

[0013] Ausgehend von einem bekannten Ansatz ergibt sich eine Modellbeschreibung, die auf einer nichtlinearen 
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Differentialgleichung basiert. Im folgenden wird eine Approximation dieser riichtlinearen Gleichung vorgestellt Im Er- 
gebnis dieser Approximation laBt sich das Systemverhatten mittels einer bilinearen Gleichung beschreiben, die die 
schnelle Losung der Beziehung im Motorsteuergerat des Kraftfahrzeugs unter Echtzeitbedingungen gestattet. Der 
gewahtte Modellansatz beinhaltet dabei die Modellierung der externen Abgasruckfuhrung und von variablen Saugsy- 
stemen und Systemen mit variabler Ventilsteuerung. Die durch diese Anordnung und durch dynamische Nachladung, 
d.h, durch Reflexionen von Druckwellen im Saugrohr hervorgerutenen Effekte, konnen ausschlieBlich durch die Wahl 
stationar bestimmbarer Parameter des Modelles sehr gut be rucks ichtigt werden. Alle Modellparameter sind einerseits 
physikalisch interpretierbar und andererseits ausschlieBlich aus stationaren Messungen zu gewinnen. 
[0014] Die meisten Algorithmen zur zeitdiskreten Losung der Differentialgleichungen lur das Model! eines Ansaug- 
systems erfordern vor allem bei geringem Druckabfall uber der Drosselklappe, d.h. bei Vollast eine sehr kleine Re- 
chenschrittweite, urn numerisch stabil zu arbeiten. Die Folge ware ein unvertretbarer Rechenaufwand bei der Bestim- 
mung der LastgroBe. Da Lasterfassungssysteme meist segmentsynchron arbeiten, d.h. fur 4-Zylindermotoren wird 
alle 180° KW ein MeBwert abgetastet, muB die Modellgleichung ebenfalls segmentsynchron gelost werden. Im nach- 
folgenden wird ein absolut stabiles Differenzen schema zur Losung von Differentialgleichungen eingesetzt, das nume- 
rische Stabilitat bei beliebiger Schrittweite garantiert. 

[0015] Dieses System zur modellgestutzten Lasterfassung ist in den bekannten Motorsteuerungssystemen, d.h. bei 
luftmassengefuhrten bzw. saugrohrdruckgefuhrten Motorsteuerungssystemen einsetzbar und im folgenden wird ein 
Korrekturalgorithmus in Form eines Model Iregelkreises formuliert, der bei auftretenden Ungenauigkeiten von Modell- 
parametern eine permanente Genauigkeitsverbesserung, d.h. einen Modellabgleich im stationaren und instationaren 
Betrieb gestattet. 

[0016] Bei der modellgestutzten Berechnung der in die Zylinder der Brennkraftmaschine einstromenden Frischluft- 
masse wird von der in Figur 1 dargestellten prinzipiellen Anordnung ausgegangen. Aus Grunden der Obersichtlichkeit 
ist dabei nur ein Zylinder der Brennkraftmaschine dargestellt. Mit dem Bezugzeichen 10 ist dabei ein Saugrohr einer 
Brennkraftmaschine bezeichnet, in dem eine Drosselklappe 11 angeordnet ist. Die Drosselklappe 11 ist miteinem, den 
Offnungsgrad dieser Drosselklappe ermittelnden Drosselklappenstellungsfuhler 14 verbunden. Stromaufwarts der 
Drosselklappe 11 ist ein Luftmassenmesser 12 und stromabwarts der Drosselklappe ist ein Saugrohrdruckfuhler 13 
im Saugrohr 10 angeordnet. Die Ausgangssignale ^d/clma/ dkw » Ps_SEN<*es Luftmassenmessers 12, des Drossel- 
klappenstellungsfuhlers 14 und des Saugrohrdrucksensors 13 sind mit Eingangen einer nicht dargestellten, an sich 
bekannten elektronischen Steuerungseinrichtung der Brennkraftmaschine verbunden. AuBerdem sind in Figur 1 noch 
schematisch ein EinlaBventil 15, ein AuslaBventil 16, sowie ein in einem Zylinder 17 beweglicher Koiben 18 dargestellt. 
[0017] zum Ruckfuhren eines Teils des Abgases in der mit einem Pfeilsymbol gekennzeichneten Richtung ist eine 
Abgasruckfuhrleitung 1 9 vorgesehen, die einen Abgastrakt 20 mit dem Saugrohr 10derart verbindet, daB Abgas strom- 
abwarts der Drosselklappe 11 in das Saugrohr 10 eingeleitet werden kann. Ein im Abgastrakt 20 angeordneter Tern- 
peratursensor 32 gibt ein der Temperatur des Abgases entsprechendes Signal an die Steuerungseinrichtung der Brenn- 
kraftmaschine ab. In die Abgasruckfuhrleitung 1 9 ist ein Abgasruckf uhrventil 21 eingeschaltet, das entweder als pneu- 
matisch arbeitendes Ventil ausgestaltet ist und uber einen elektropneumatischen Wandler mittels eines Tastverhalt- 
nisses angesteuert wird oder als elektrisch betatigtes Ventil realisiert ist und durch entsprechende Ansteuersignale 
vom elektronischen Steuergerat der Brennkraftmaschine den Offnungsquerschnitt der Abgasruckfuhrleitung 19 ver- 
andert Dem Abgasruckfuhrventil 21 ist ein Stellungsf uhler 22 zugeordnet der ein dem Offnungsgrad des Abgasruck- 
fuhrventils 21 entsprechendes Signal AGRW an die elektronische Steuerungseinrichtung abgibt. 
[0018] AuBerdem sind in Figur 1 ausgewahlte GroBen bzw. Parameter des Saugsystems eingezeichnet. Dabei be- 
deutet das Dachsymbol "a" uber einer GroBe, daB es sich urn eine ModellgroBe handelt, wahrend GroBen ohne Dach- 
symbol MeBgroBen reprasentieren. Im einzeinen bedeuten: 

[0019] P y Umgebungsdruck, P s Saugrohrdruck, T s Temperatur der Luft im Saugrohr, V s Saugrohrvolumen, P A Ab- 
gasgegendruck, T A Abgastemperatur, P RG Restgaspartialdruck im Saugrohr, P FG Frischgaspartialdruck im Saugrohr. 
[0020] GroBen mit einem Punktsymbol kennzeichnen die erste zeitliche Ableitung der entsprechenden GroBen. Mit 
m DK ist somit der Luftmassenstrom an der Drosselklappe, mit m RG der Restgasmassenstrom am Abgasruckfuhrventil, 
mit m Zyf FG der Frischgasmassenstrom in den Zylinder und mit m 2yl RG der Restgasmassenstrom und mit m Zyf der 
Gesamtmassenstrom in den Zylinder bezeichnet. 

[0021] Sind der Abgasgegendruck, der Offnungsgrad des Abgasruckf uhrventils und die Abgastemperatur bekannt, 
entweder durch eine betriebspunktabhangige Parameterisierung mittels Kennfelder oder durch Auswerten der Signale 
der entsprechenden Sensoren, so kann der Massenstrom durch das Abgasruckfuhrventil berechnet werden. 
[0022] Unter der Voraussetzung konstanter Temperatur der Luft im Saugrohr laBt sich die Differentialgleichung aus 
der Zustandsgleichung idealer Gase zur Berechnung des Absolutdrucks im Saugrohr aufstellen (vgl. Fig.1) : 
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(mDK+mRG-mxyt) 
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[0023] Dabei gilt: 

h zeitliche Ableitung der ModellgroBe Saugrohrdruck 
allgemeine Gaskonstante 
Temperatur der Luft im Saugrohr 
Saugrohrvolumen 
iox : Modelluftmassenstrom an der Drosselklappe 
lt K : Modellrestgasmassenstrom am Abgasruckfuhrventil 
"m» : ModellgroBe des Gesamtmassenstroms in den Zylinder 

[0024] Die durch (1.0) beschriebenen Verhaltnisse sind auf Mehrzylinder-Brennkraftmaschinen mit Schwingrohr- 
(Schaltsaugrohr-) und/oder Resonanzsaugsysteme ohne strukturelle Anderungen anwendbar. 
[0025] Fur Systeme mit Multi-Point-Einspritzungen, bei denen die Kraftstoffzumessung durch mehrere Einspritzven- 
tile erfolgt, gibt die Gleichung (1.0) die Verhaltnisse genauer wieder afs dies bei Single-Point-Einspritzungen, d.h. bei 
Einspritzungen, bei denen der Kraftstoff mittels eines einzigen Kraftstoff einspritzventiles zugemessen wird, der Fall 
ist. Bei erstgenannter Art der Kraftstoffzumessung ist nahezu das gesamte Ansaugsystem mit Luft gefullt. Lediglich in 
einem kleinen Bereich vor den EinlaGventilen befindet sich ein Kraftstoff-Luftgemisch. Im Gegensatzdazu ist bei Single- 
Point-Einspritzsystemen das gesamte Saugrohr von der Drosselklappe bis zum EinlaBventil mit Kraftstoff- Luft-Ge- 
misch gefullt: da das Einspritzventil vor der Drosselklappe angeordnet ist. In diesem Fall stellt die Annahme eines 
idealen Gases eine starkere Naherung dar, als dies bei der Multi-Point-Einspritzung der Fall ist. Bei Single-Point- 
Einspritzung erfolgt die Kraftstoffzumessung entsprechendiu* , bei Multi-Point-Ein-spritzungen entsprechendm*. 
[0026] Im folgenden wird die Berechnung der Massenstromei,^ undi^, naher beschrieben. 
[0027] Die ModellgroBe des Luftmassenstromes an der Drosselklappe^ wird durch die DurchfluBgleichung idealer 
Gase durch Drosselstellen beschrieben. Der Luftmassenstronrw wird demnach durch die Beziehung 
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mit der DurchfluBfunktion 
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fur Oberkritische Druckverhaltnisse bzw. 
V = const, fur kritische Druckverhaltnisse 
bestimmt. 
so [0028] Dabei bedeuten: 



ModellgroBe des Luftmassenstromes an der Drosselklappe 


reduzierter Strdmungsquerschnitt 

k: 

Adiabatenexponent 

R L : 

allgemeine Gaskonstante 

Is- 

Temperatur der Luft im Saugrohr 

Zu- 

ModellgroBe des Umgebungsdruckes 

Ps- 

ModellgroBe des Saugrohrdruckes 
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DurchfluBfunktion 


10 


[0029] An der Drosselstelle, d.h. an der Drosselklappe 11 auftretende Strom ungsverluste werden durch den redu- 
zierten Qu rschnitt an der Drosselklappe A^q^ berucksichtigt. Aus stationaren Messungen kann bei bekannten DrGk- 
ken vor und hinter der Drosselstelle und bekanntem Massenstrom durch die Drosselstelle eine Zuordnung zwlschen 
dem vom Drosselklappenstellungsfuhler 14 ermittelten Drosselklappenwinkel DKW und dem entsprechendem redu- 
zierten Querschnitt A TOK ^angegeben werden. 

[0030] Wlrd der Luftmassenstronw an der Drosselklappe durch die Beziehung (2.0) beschrieben, so entsteht ein 
komplizierter Algorithmus zur numerisch richtigen Losung der Differentialgleichung (1.0) mit einem absolut stabilen 
Differenzenschema. Zur Reduktion des Rechenaufwandes wird die DurchfluBfunktion \|/ durch einen Polygonzug ap- 
proximiert. 

[0031] Figur 2 zeigt den Verlauf der DurchfluBfunktion y und das darauf angewandte Approximationsprinzip. Inner- 
halb eines Abschnittes i,j (i.j = 1 ...k) wird die DurchfluBfunktion \jr durch eine Gerade dargestellt. Mit einer vertretbaren 
Anzahl von Geradenabschnitten kann damit eine gute Approximation der DurchfluBfunktion durch 


15 


20 


A 


(1. . .k) 
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erreicht werden. 

[0032] In dieser Form beschreibt m { die Steigung und t\ { das Absolutglied (Offset) des jeweiligen Geradenabschnittes. 
In den einzelnen Abschnitten der Unearlslerung gellen jeweils die Werte m j} n t Die Werte fur die Steigung und fur das 
Absolutglied werden in Tabellen als Funktion des Verhaltnisse^ Saugrohrdruck zu Umgebungsdruck-^ abgelegt. 
[0033]a Auf der Abszisse von Figur 2 ist das Druckverhaltriis-*^ bzw. das Verhaltnis Saugrohrdruck ziFAbgasgegen- 
druck-*^ und auf der Ordinate der Funktionswert (0 - 0.3) der CMjrchfluBfunktion y aufgetragen. 
[0034] *Durch einen solchen Ansatz kann die Gleichung (2.0) zur Berechnung des Luftmassenstromes an der Dros- 
selklappe^ durch die Beziehung 
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fur i = (1 ...k) 
approximiert werden. 

[0035] Die Berechnung des Restgasmassenstroms am AbgasrGckfGhrventil erfolgt ebenfalls mit der DurchfluBglei- 
chung idealer Gase fur Drosselstellen. Am AbgasruckfGhrventil herrscht jedoch das Druckgefallezwischen Abgasge- 
gendruck P A und Saugrohrdruck p s , sowie das Abgastemperaturniveau T A . Fur die Approximation der DurchfluBfunk- 
tion y sind in diesem Fall die Werte mj und nj zu berucksichtigen, da sich der Umgebungsdruck P u vom Abgasgegen- 
druck P A unterscheidet (Figur 2). Auftretende Strdmungsverluste werden durch den reduzierten Stromungsquerschnitt 
an der Drosselstelle AbgasruckfGhrventil A RRG beschrieben. 


so 
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[0036] Die in den jeweiligen Zylinder der Brennkraftmaschine einstromende Luftmasse laBt s\ch analytisch nut 
schwer berechnen, da sie stark vom Ladungswechsel abhangt. Die Fullung der Zylinder wird weitestgehend durch den 
Saugrohrdruck, die Drehzahl und durch die Ventilsteuerungszeiten bestimmt. 

[0037] Zur moglichst genauen Berechnung des Massenstroms in den jeweiligen Zylinder^ ist deshalb einerseits 
die Beschreibung der Verhaltnisse im Ansaugtrakt der Brennkraftmaschine mittels partieller Differentialgleichungen 
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und andererserts die Berechnung des Massenstromes am EinlaBventi! nach der DurchffuBgfeichung als erforderfiche 
Randbedingung notwendig. Erst dieser komplizierte Ansatz gestattet die Berucksichtigung dynamischer.Nachladeef-. 
fekle, die von der Drehzahl, der Saugrohrgeometrie, der Zylinderzahl sowie den Ventilsteuerzeiten maBgeblich beein- 
fluBt werden. 

[0038] Da eine Berechnung nach oben genanntem Ansatz in der elektronischen Steuerungseinrichtung der Brenn- 
kraftmaschine nicht realisierbar ist, geht eine mogliche Naherung von einem. einfachen Zusammenhang zwischen 
Saugrohrdruck P s und Zyindermassenstrornii*, aus. Fur einen weiten Bereich der sinnvollen Ventilsteuerzeiten kann 
dafur in guter Naherung von einem linearen Ansatz der Form 


(4.0) 


ausgegangen werden. 

[0039] Die Steigung y 1 und das Absolutglied Yo der Beziehung (4.0) sind dabei, unter Berucksichtigung aller wesent- 
lichen EinfluBfaktoren Funktionen der Drehzahl, der Saugrohrgeometrie, der Zylinderzahl, der Ventilsteuerzeiten sowie 
der TemperaXur der Luft im Saugrohr T s . Die Abhangigkeit der Werte von Yf und y 0 von den EinfJuBgrdBen Drehzahl, 
Saugrohrgeometrie, Zylinderzahl und den Ventilsteuerzeiten und Ventilerhebungskurven kann dabei uber stationare 
Messungen ermittelt werden. Uber diese Wertebestimmung wird ebenfalls der EinluB von Schwingrohr- und/oder Re- 
sonanzsaugsystemen auf die von der Brennkraftmaschine angesaugte Luftmasse gut wiedergegeben. Die Werte von 
und yo sind in Kennfeldern der elektronischen Motorsteuerungseinrichtung abhangig von der Drehzahl n abgelegt. 
[0040] Als bestimmende GrdBe zur Ermittlung der Motorlast wird der Saugrohrdruck P s ausgewahlt. Mrt Hilfe der 
Modell-Differentialgleichung sod diese GrdBe mdglichst exakt und schnell geschatzt werden. Die Schatzung von P s 
erfordert die Losung der Gleichung (1.0). 

[0041] Mit den Gleichungen (2.3), (3.0) und (4.0) kann die Formel (1.0) durch die Beziehung 
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approximiert werden. 

[0042] Zur Losung der Gleichung (5.0) wird diese Beziehung in eine geeignete Differenzengleichung ubergefuhrt. 
[0043] Als Kriterium zur Auswahl des geeigneten Differenzenschemas konnen die tolgenden prinzipiellen Anforde- 
rungen an die Losungseigenschaften der zur bildenden Differenzengleichung formuliert werden: 

1. Das Differenzenschema muB auch unter extremen dynamischen Anforderungen konservativ sein, d.h. die Lo- 
sung der Differenzengleichnung muB der Losung der Differentialgleichung entsprechen, 

2. die numerische Stabilitat muB fur Abtastzeiten, die den maximal moglichen Segmentzeiten entsprechen, im 
gesamten Arbeitsbereich (Dynamikbereich) das Saugrohrdruckes garantiert sein. 

[0044] Forderung 1 ist durch einen impliziten Rechenalgorithmus erfuHbar. Aufgrund der Approximation der nichtli- 
nearen Differentialgleichung (1.0) durch eine bilineare Gleichung (5.0) ist das entstehende implizite Losungsschema 
ohne Einsatz iterativer Verfahren losbar, da die Differenzengleichung in eine explizite Form Qberf uhrt werden kann. 
[0045] Die zweite Forderung ist aufgrund der Konditionierung der Differentialgleichung (1.0) und deren Approxima- 
tion (5.0) nur durch e\ne Rechenvorschrift zur Bildung der Differenzengleichung erfullbar, die absolut stab'il arbeitet. 
Diese Verfahren werden auch als A-stabile Verfahren bezeichnet. Kennzeichnend fur diese A-Stabilitat ist die Eigen- 
schaft des Algorithmus, bei einem stabilen Ausgangsproblem fur beliebige Werte der Abtastzeit, d.h. Segmentzeit T AB 
numerisch stabil zu sein. Eine mogliche Rechenvorschrift zur numerischen Losung von Differentialgleichungen, die 
beiden Forderungen gerecht wird, ist die Trapezregel. 

[0046] Die durch Anwendung der Trapezregel entstehende Differenzengleichung lautet im vorliegenden Fall 
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(5.1) 


furN = (L..oo) 

[0047] bezeichnet die Zeit zwischen 2 Abtastschritten N. 

[0048] Wird die D iff e rent ialgleichung (5.0) mit Hilfe der Trapezregel (5.1 ) gelost, so ergibt sich fur den Saugrohrdruck 
10 zum Abtastschritt N folgende Beziehung. 


75 


20 


25 


P S [N] = 


V 


2X 1 


2X 1 


2% 1 


fur N « (l...oo) und i,j = (l...k) 


(5.2) 


[N] bedeutet dabei das aktuelle Segment bzw. den aktuellen Rechenschritt, [N+1] das nachstfolgende Segment bzw. 
30 den nachstfolgenden Rechenschritt. 

[0049] 1st der Gesamtdruck im Saugrohr bekannt, so kann die Anderung des Restgaspartialdrucks P RG im Saugrohr 
nach foigender Beziehung berechnet werden: 


35 


40 


45 


R -T • 
P RC — (mRc-mzyi_*c) 


mit der Gleichung (3.0) 


als Restgasmassenstrom durch das Abgasruckfuhrventii und mit 


(6.0) 


50 


55 


(7.0) 


als Restgasmassenstrom in den Zylinder durch externe Abgas rue kfuh rung. 

[0050] WendetmandieTrapezregel (5.1) zur Losung der Differentialgleichungfur den Restgaspartialdruckim Saug- 
rohr (6.0) an, so ergibt sich die Beziehung 
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v s V 


2 k ) 


(8.1) 


[0051] zur Berechnung des Restgaspartialdruckes zum Zeitpunkt N. Weiterhin gilt die Beziehung: 

P S = Pf B +Pr B (90) 


A 

[0052] Der Frischgaspartialdruck P FG ist dann: 

Pfg-Ps-P^ < 9 - 1 ) 

[0053] Damit sind der Frischgas - und Restgaspartialdruck irruSaugrohr im stationaren und instationaren Motorbe- 
trieb bekannt. Mit dem Frischgaspartialdruck kann der Frischluftmassenstrom in die jeweiligen Zylinder berechnet 
werden: 


™ZYL FG=1\' P FG+t* (10.0) 


[0054] Durch einen einfachen Integrationsalgorithmus erhalt man die wahrend eines Ansaugtaktes von der Brenn- 
kraftmaschine angesaugte Luftmasse. 


m 


l ZYL _FG 


= "f - • (™zn _fg[X] + m ZYL _ FC [ N - 1]) (10.1) 


[0055] Durch systembedingte Totzeiten bei der Kraftstoffdosierung und einer notwendigen Kraftstoffvorl age rung ist 
eine Vorhersage der von der Brenn kraftmaschine angesaugten Luftmasse mit einem wahlbaren Vorhersagehorizont 
wunschenswert, da nur so im instationaren Betrieb eine genaue Einhaltung eines geforderten Kraftstoff-Luft-Verhalt- 
nisses moglich ist. 

[0056] Bei dem beschriebenen Verfahren liegen die ModellgrdGen der zeitlichen Anderung der Partialdrucke im An- 
saugsystem in analytischer Form vor. Durch H-fache Anwendung der Trapezregel kann der Frischgaspartialdruck durch 
die Beziehung 


(10.2) 

um H Segmente vorhergesagt werden. 

[0057] Geht man von uber die Vorhersagezeif konstanten Werten und Yo aus, so kann die vorhergesagte Frisch- 
gasmasse im Zylinder durch die Gleichung 
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(10.3) 

to bestimmt werden. 

[0058] Wahlt man fur den Vorhersagehorizoht (Pradiktionshorizont) H Werte in der GroBenordnung von 1....3, so 
konnen mit der Forme! (1 0.3) Vorhersageergebnisse fur die von der Brennkraftmaschine angesaugten Frischgasmasse 
mit hoher Genauigkeit erhalten werden. 

[0059] Im folgenden wird das Prinzip des Modellabgleichs fur luftmassen- und saugrohrdruckgefOhrte Motorsteue- 
*5 rungssysteme beschrieben. 

[0060] . Bedingt durch den Einsatz von Motoren mit variabler Ventilsteuerung und/oder veranderlicher Saugrohrgeo- 
metrie, durch Fertigungstoleranzen und Alterungserscheinungen, sowie durch Temperatureinflusse sind die Werte von 
y n und Y 0 mit einer gewissen Unsicherheit behaftet. Die Parameter der Gleichung zur Bestimmung des Massehstromes 
in den Zylindem sind, wie oben beschrieben, Funktionen vielfaltiger EinfluBgroBen, von denen nur die wichtigsten 
20 erfaBt werden konnen. 

[0061] Bei der Berechnung des Massenstromes an der Drosselklappe wirken sich MeBfehler bei der Erfassung des 
Drosselklappenwinkels und Approximationsfehler bei der Polygonzugapproximation der DurchfluBfunktion \j/ auf die 
ModellgroBen aus. Besonders bei kleinen Drosselklappenwinkeln ist die Systemempfindlichkeit gegenuber erstge- 
nannten Fehlern besonders hoch. Daraus ergibt sich : daB kleine Anderungen der Drosselklappenstellung einen gra- 
25 vierenden EinluB auf Massenstrom bzw. Saugrohrdruck haben. Urn die Wirkung dieser Einflusse zu reduzieren, wird 
im folgenden ein Verfahren vorgeschlagen, das es gestattet, bestimmte GroBen, die EinfluB auf die Modellrechnung 
haben; so zu korrigieren, daB eine genauigkeitsverbessernde Mode I Ian pass ung fur stationaren und instationaren Mo- 
torbetrieb durchgefuhrt werden kann. 

[0062] Die Anpassung der wesentlichen Parameter des Modells zur Bestimmung der LastgroBe der Brennkraftma- 
30 schine erfolgt durch die Korrektur des aus dem gemessenen Drosselklappenwinkel bestimmten reduzierten Stro- 
mungsquerschnittes A RDK durch die KorrekturgroBe AA RDK : 

A RDK_KORR ~ A RDK + ^ A RDK ( 1 1 °) 

35 

[0063] Je nachdem welcher Lasterfassungssensor verwendet wird, erfolgt die Berechnung der KorrekturgroBe 
AA RDK in einem Modellregelkreis. Bei einer Luftmassenmessung erfolgt die Berechnung der Regelabweichung zwi- 
schen dem gemessenen Luftmassenstrom an der Drosselklappe und dem Model luftmassenstrom an der Drosselklap- 
pe. 

40 [0064] Bei einer Saugrohrdruckmessung erfolgt die Berechnung der Regelabweichung zwischen dem Modellsaug- 
rohrdruck und dem gemessenen Saugrohrdruck. 

[0065] Fur luftmassengefuhrte Motorsteuerungssysteme ist also der mittels des Luftmassenmessers an der Dros- 
selklappe gemessene Luftmassenstrom rnoK_LMN d ' e FuhrungsgroBe dieses Regelkreises, wahrend fur saugrohr- 
druckgefOhrte Systeme der mittels des Saugrohrdrucksensors gemessene Saugrohrdruck P s s^als FuhrungsgroBe 
45 genutzt wird. Uber eine Folgeregelung wird der Wert von AA RED dann so bestimmt, daB die Regelabweichung zwischen 
FuhrungsgroBe und der ensprechenden RegelgroBe minimiert wird. 

[0066] Um auch im dynamischen Betrieb Genauigkeitsverbesserungen mit der genannten Methode zu erreichen, 
muB die MeBwerterfassung der FuhrungsgroBen moglichst exakt nachgebildet werden. In den meisten Fallen sind 
dabei das dynamische Verhalten des Sensors, d.h. entweder des Luftmassenmessers oder des Saugrohrdrucksensors 
50 und eine nachfolgend durchgefuhrte Mittelwertbildung zu berucksichtigen. 

[0067] Das dynamische Verhalten des jeweiligen Sensors kann in erster Naherung als ein System erster Ordnung 
mit eventuell arbeitspunktabhangiger Verzogerungszeit Tj modelliert werden. 
■ [0068] Im Falle einer saugrohrdruckgefuhrten Motorensteuerung lautet eine mdgliche Gleichung zur Beschreibuhg 

des Sensorverhaltens 

55 
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hjBiW* h.snlv-ih T ' +P s [N-i]\ i 


(12.0) 


wahrend fQr ein luftmassengefuhrtes System die Beziehung 


( 


(13.0) 


angegeben werden kann. 7 n bezeichnet die Verzogerungszeit des jeweiligen Sensors. 

[0069] Die im Steuergerat durchgefuhrte Mittelwertbildung der MeGwerte uber ein Segment kann durch die Bezie- 
hungen 

oder 


**DK_U*4 ,MTT [ N] = ^ ' {™DK-IMM [ N ~\] + *l DK _ ^ [ AT]) ( 1 3 . 1 ) 


modelliert werden. 

[0070] Falls beide Lasterf assungssensoren vorhanden sind, kann zum Abgleich des Luftmassenstroms an der Dros- 
selklappe mit dem gemessenen Luftmassenstrom der reduzierte Stromungsquerschnitt an der Drosselklappe bzw. der 
Umgebungsdruck verwendet werden. Der Saugrohrdrucksensor kann im Abgasruckfuhrungsbetrieb zur Bestimmung 
des reduzierten Stromungsquerschnittes am Abgasruckfuhrungsventil bzw. des Abgasgegendrucks verwendet wer- 
den. 

[0071] Treten Differenzen zwischen dem Modellsaugrohrdruck und dem gemessenen Saugrohrdruck auf, so werden 
sie durch die Anpassung des Restgasmassenstroms uber die Korrektur A RRG durch AA RRG ausgeregelt (Figur 3). 
Uberschreitet die StellgroBe AA RRG eine bestimmte Schwelle in positiver Richtung so erfolgt eine Erhohung der Mo- 
dellgroGe Abgasgegendruck, wird eine negative applizierbare Schwelle unterschritten : so erfolgt eine Erniedrigung der 
ModellgroGe Abgasgegendruck. 

[0072] Bei geschlossenem Abgasruckfuhrventil wird mit dem Saugrohrdrucksensor der Zusammenhang zwischen 
Saugrohrdruck und Luftmasse mit Hilfe des Parameters y 1 korrigiert (Figur 4). Die KorrekturgroGe wird drehzahlabhan- 
gig bestimmt und im nichtfluchtigen Speicher des Steuergerats abgelegt. Der Korrekturwert wird auf einen Maximalwert 
begrenzt. 

[0073] im folgenden wird ein Modellabgleich fur Motorsteuerungssysteme beschrieben, wenn beide Lasterf assungs- 
sensoren (Luftmassenmesser und Saugrohrdrucksensor) vorhanden sind. Fur dieses System konnen die in den Figu- 
ren 3 und 4 dargestellten Modellstrukturen angegeben werden. 

[0074] Der Drosselklappenstellungsfuhler 14 (Figur 1) liefert ein dem Offnungsgrad der Drosselklappe 11 entspre- 
chendes Signal, z.B. einen Drosselklappenoffnungswinkel DKW. In einem ersten Kennfeld KF1 der elektronischen 
Motorsteuerungseinrichtung sind zu verschiedenen Werten dieses Drosselklappenoffnungswinkels zugehorige Werte 
fur den reduzierten Querschnitt der Drosselklappe A RDK abgespeichert. Das Teilsystem 'SaugrohrmodeH" in den Fi- 
guren 3 und 4 reprasentiert das durch Gleichungen (5.2) und (8.1) beschriebene Verhalten. 

[0075] FuhrungsgroBen der Modellregelkreise sind die Me3werte der uber ein Segment gemittelten GroGen des 
Lutmassenstromes an der Drosselklappe^ ^ und des Saugrohrdruckes P s _sen- Werden als Regler in den Modell- 
kreisen Pl-Reglereingesetzt, soistdie bleibende Regelabweichung null, d.h. im stationaren Fall sind die ModellgroGen 
und die entsprechenden MeGgroGen identisch. 

[0076] Die Pulsationserscheinungen des Luftmassenstromes an der Drosselklappe, die vor allem bei 4-Zylindermo- 
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toren zu beobachten sind, fQhren bei betragsbildenden Luftmassenmessern zu erheblichen positiven MeBfehlem und 
somit zu einer stark fehlerbehafteten FuhrungsgroBe. Durch eine Abschaltung des entsprechenden Reglers, d.h. einer 
Verkleinerung der Reglerparameter kann zum gesteuerten modellgestutzten Betrieb Qber; gegangen werden. Berei- 
che, in denen die genannten Pulsationen auftreten, konnen somit mit dem selben Verfahren unter Berucksichtigung 

s dynamischer Zusammenhange behandelt werden, wie diejenigen Bereiche, in denen eine nahezu ungestorte Fuh- 
rungsgroBe vorliegt. Im Gegensatz zu Verfahren, die relevante MeBwerte nur in stationaren Betriebspunkten beruck- 
sichtigen, bleibt das beschri bene System nahezu uneingeschrankt arbeitsfahig. Bei Ausfall des Signals des Luftmas- 
senmessers oder des Drosselklappenstellungsfuhlers ist das vorgestellte System in der Lage, ein entsprechendes 
Ersatzsignal zu bilden. Bei Ausfall der FuhrungsgroBe muB der gesteuerte Betrieb reallsiert werden, wahrend im an- 

10 deren Fall der geregelte Betrieb die kaum beeintrachtigte Funktionsfahigkeit des Systems garantiert. 

[0077] Der Block ■Saugrohnmodell" reprasentiert die Verhaltnisse wie sie anhand der Gleichungen (5.2) und (8.1) 
beschrieben sind und hat demzufolge als AusgangsgrdBen die ModellgroBen P s und th M , sowie den Wert fur den 
Frischgasmassenstrom'^w.re in den Zylinder der Brennkraftmaschine, der als Grundlage zur Bestimmung der Ba- 
siseinspritzzeit herangezogen wird. Nach der Modellierung des Sensorubertragungsvemaltens des Luftmassenmes- 

15 sers und des Saugrohrdrucksensors (Gleichungen 1 2.0 bzw. 1 3.0) werden die ModellgroBen P s und^ einer Mittelung 
gemaB Gleichung (12.1) bzw. (13.1) unterzogen, so daB die gemittelten GroBen \s£A/_M/r»Vwf.i«' unddervom 
Saugrohrdrucksensor gemessene Wert Ps_sen dzw der v ° m Luftmassenmesser gemessene Luftmassenstrom 
m DK jjMf jeweils einer Vergleichsstelle 23, 24 zugefuhrt werden konnen. 

[0078] Die Differenz zwischen der MeBgroBe und der ModellgroBe des gemittelten Massenstromes an der Drossel- 
20 klappe, die an der Vergleichstelle 24 gebildet wird, dlent in Figur 3 und Figur 4 als Grundlage zur Berechnung des 
Korrekturwertes AA RDK '\m Regler 27. Die ModellgroBe des reduzierten Querschnittes wird entsprechend der Gleichung 
(11.0) gebildet und gestattet die Korrektur der ModellgroBe des Massenstromes an der Drosselklappe im stationaren 
und instationarem Motorbetrieb. 

[0079] In Figur 3 wird die Differenz zwischen der gemittelten MeBgroBe und der entsprechenden ModellgroBe des' 
25 Saugrohrdruckszur Berechnung eines Korrekturwertes AA RRG zur Beeinflussungdes Restgasmassenstromes an dem 
AGR-Ventil verwendet. Uberschreitet die KorrekturgroBe AA RRG eine in einer Schwellwertstule 29 festgelegte Schwel- 
le, so erfolgt die Korrektur des Abgasgegendruckes P A in positiver bzw. negativer Richtung entsprechend des Vorzei- 
chens von AA RRG . 

[0080] In Figur 4 wird die Differenz zwischen der gemittelten MeBgroBe und der entsprechenden ModellgroBe des 
30 Saugrohrdruckes zur Korrektur des Massenstromes in den Zylinder durch eine Beeinflussung von y A verwendet. 

[0081] Das in Figur 3 vorgestellte Verfahren bietet den Vorzug, daB auch bei Parameterfehlern des Modells zur 
LastgroBenschatzung glatte Verlaufe der LastgroBe zu erwarten sind, da die Differenzenordnung zwischen dem Stell- 
eingriff und der LastgroBe eins betragt. Allerdings ist eine Schatzung des Umgebungsdruckes in Abhangtgkeit des 
Wertes von A^A^^bzw. des Druckverl^altnisses erforderlich. Das in Figur 4 vorgestellte Verfahren bietet den Vorzug, 
35 auch bei groBen Druckverhaltnissen den Massenstrom in den Zylinder in weiten Bereichen andem zu konnen. 
Veranderungen der DurchfluSbeiwerte fefer EinlaBventile sind mit diesem Verfahren gut korrigierbar. 


Patentanspruche 

40 

1. Verfahren zum Bestimmen der einstromenden Frischluftmasse in den Zylinder einer Brennkraftmaschine mit 

einem Saugrohr (10) und eine darin angeordnete Drosselklappe (11), sowie eine den Offnungsgrad (DKW) 
der Drosselklappe (19) erfassende Sensoreinrichtung (14), 
45 - einer Abgasruckfuhreinrichtung (19) mit einem Abgasruckf Qhrventil (21 ), sowie ein den Offnungsgrad (AGRW) 

des Abgasruckfuhrventils (21) erfassenden Sensoreinrichtung (22), 

einem, ein Lastsignal (rnoK_LMM> p s_sEN) der Brennkraftmaschine erzeugenden Sensor (12; 13) 

einer elektronischen Steuerungseinrichtung, die auf der Grundlage des gemessenen Lastsignals {m DK LMh/ j, 

p s_sen) und der Drehzahl (n) der Brennkraftmaschine eine Grundeinspritzzeit berechnet, 

so 

dadurch gekennzeichnet, daB 

. . - die Verhaltnisse im Saugrohr (10) mittels eines Saugrohrfultungsmodells (Gleichungen 5.2,8.1) nachgebildet 
werden, wobei als EingangsgroBen des Modells zumindest der Offnungsgrad (DKW) der Drosselklappe (11), 
55 der Offnungsgrad (AGRW) des Abgasruckfuhrungsventils (21 ), der Umgebungsdruck (P u ), die Abgastempe- 

ratur (T A ), die Temperatur im Saugrohr (T s ) und die Ventilstellung reprasentierende Parameter herangezogen 
werden, 

jeweils eine ModellgroBe fur den Luftmassenstrom U M ) an der Drosselklappe (11) und fur den Restgasmas- 
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senstrom (*«) am AbgasruckfOhrventil (21) gebildet wifd (Gleichung 2.3; 3.0) 

jeweils eine ModellgroBe fur den Luftmassenstrom in den Zylinder bzw. fur den Restgasmassenstrom in 
. den Zylinder (^.« ) als Fun kt ion des Saugrohrdrucks (P s ) bzw. des Restgaspartialdruckes (Pro) gebildet wird 
(Gleichung 4.0; 7.0), 

aus den ModellgrdBen Luftmassenstrom an der Drosselklappe (£«), Restgasmassenstrom am Abgasruck- 
fOhrventil (i K ) und Luftmassenstrom in den Zylinder (^) als bestimmende GrdGe zur Ermittlung der tatsach- 
lichen Last der Brennkraftmaschine der Saugrohrdruck (P s ) berechnet wird (Gleichung 5.2), 
aus den Model IgroBen Restgasmassenstrom am AbgasruckfOhrventil (^J und Restgasmassenstrom in den 
Zylinder (<U.«) der Restgaspartialdruck {P RG ) im Saugrohr berechnet wird (Gleichung 8.1) 
aus dem Saugrohrdruck (P s ) und dem Restgaspartialdruck (Pr G ) der Frischgaspartialdruck (P FG ) bestimmt 
und 

aus dem Zusammenhang zwischen Frischgaspartialdruck (P FG ) und Frisch luftmassenstrom in den Zylinder 
(i*,.**) durch Integration die in den Zylinder (17) einstrdmende Luftmasse (nr*zyLFG) ©rhalten wird. 

Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daB die ModellgroBen fur den Luftmassenstrom (£«) an der Drosselklappe (11 ) und fur den Restgasmassenstrom 
(*«) am AbgasruckfOhrventil (21) mit Hilfe der DurchfluBgleichung idealer Gase durch Drosselstellen beschrieben 
werden. 

Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daG an den Drosselstellen (11 ,21 ) auftretende Sromungsverluste jeweils durch die Verwendung reduzierter Strd- 
mungsquerschnitte (A RDK , A RRG ) an den Drosselstellen be rucks ichtigt werden. 

Verfahren nach Anspruch 3, 
dadurch gekennzeichnet, 

daG die Werte f Or die reduzierten Querschnitte (A rdk ,Arr G ) aus stationaren Messungen abhangig vom Offnungs- 
grad (AGRW) des AbgasrOckfOhrventils am MotorprOfstand festgelegt und in Kennfeldem (KF1; KF4) eines Spei- 
chers der Steuerungseinrichtung abgelegt sind , 

Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daG bei der Darstellung der ModellgroBen fOr den Luftmassenstrom an der Drosselklappe (11) und der Mo- 
dellgroBe fur den Restgasmassenstrom (l K ) am AbgasruckfOhrventil (21) eine in den DurchfluBgleichungen vor- 
handene DurchfluBfunktion (y) in einzelne Abschnitte (i,j = 1...k) unterteilt wird und diese Abschnitte durch Gera- 
denabschnitte angenahert werden, wobei fur die Werte der Steigung (nrij j) und f Or da^ Absolutglied (n L |) der jewei- 
ligen Geradenabschnitte jeweils als Funktion des^erhaltnisses von Saugrohrdruck (P s ) zu Umgebungsdruck (P u ) 
bzw. Saugrohrdruck (P s ) zu Abgasgegendruck [P A ) bestimmt und in einem Kennfeld eines Speichers der Steue- 
rungseinrichtung abgelegt sind. 

Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daB die ModellgroBe fOr den Luftmassenstrom in den Zylinder als eine lineare Funktion des Saugrohrdrucks 
(P s ) beschrieben wird (Gleichung 4.0) und die Steigung (y^ und das Absolutglied (y 0 ) der linearen Funktion ab- 
hangig von mindestens einem der Parameter, Drehzahl der Brennkraftmaschine, Zylinderzahl, Saugrohrgeome- 
trie, Temperatur der Luft (T s ) im Ansaugrohr (10) und Ventilsteuerzeiten festgelegt sind. 

Verfahren nach Anspruch 1 , 
dadurch gekennzeichnet, 

daB die ModellgroBe fur den Restgasmassenstrom in den Zylinder (m WJW ) als eine lineare Funktion des Restgas- 
partialdrukkes (P RG ) beschrieben wird (Gleichung 7.0) und die Steigung (y^ der linearen Funktion abhangig von 
mindestens einem der Parameter, Drehzahl der Brennkraftmaschine, Zylinderzahl, Saugrohrgeometrie, Tempe: 
ratur der Luft (T s ) im Ansaugrohr (10) und Ventilsteuerzeiten festgelegt sind. 

Verfahren nach Anspruch 6 oder 7, 
dadurch gekennzeichn t, 

daB die Parameter durch stationare Messungen am MotorprOfstand ermittelt werden und in Kennfeldem (KF2,KF3) 
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eines Speichers der Steuerungseinrichtung abgelegt sind. : 


9. Verfahren nach Anspruch 1, 
dadurch g kennzeichnet, 

daB das vom Lastsensor (1 2; 1 3) gemessene Lastsignal (m DK P s _sen) zur Korrektur und damit zum Abgleich 
der ModellgrdBen Luftmassenstrom an der Drosselklappe (jL~) und Saugrohrdruck (P^ in einem Modellregelkreis 
herangezogen wird, wobei das Lastsignal (rn DK LMA/ f, Ps_sen) a,s FuhrungsgroBe des Regelkr eises dient 

10. Verfahren nach Anspruch 9, 
dadurch gekennzeichnet, 

daB der Abgleich im stationaren und/oder instationaren Betrieb der Brennkraftmaschine durchgefuhrt wird und 
dabei das Obertragungsverhaiten (Gleichung 12.0; 13.0) des Lastsensors (12; 13) berucksichtigt wird. 

11. Verfahren nach Anspruch 10, 
dadurch gekennzeichnet, 

daB dais Ubertragungsverhatten des Lastsensors (12,13) durch ein System erster Ordnung mit arbeitspunktab- 
hangiger Verzogerungszeit (T1 ) modelliert wird. 

12. Verfahren nach Anspruch 10, 
dadurch gekennzeichnet, 

daB jedem gemessenen Wert des Drosselklappenoffnungsgrades (DKW) ein Wert eines reduzierten Querschnitts 
der Drosselklappe (A RDK ) zugeordnet ist und der Abgleich der ModellgroBe Luftmassenstrom an der Drosselklappe 
(l„ ) durch Korrektur des reduzierten Querschnitts (A RDK ) durch eine KorrekturgroBe (AA RDK ) derart erfolgt, daB 
die Regelabweichung zwischen FuhrungsgroBe (rn DK LMM ) und der ModellgroBe (L) minimiert wird. 

13. Verfahren nach Anspruch 10, 
dadurch gekennzeichnet, 

daB jedem gemessenen Wert des Offnungsgrades (AGRW) des Abgasruckfuhrventils (21 ) Wert eines reduzierten 
Querschnitts der des Abgasruckfuhrventils {A RRG ) zugeordnet ist und der Abgleich der ModellgroBe Saugrohrdruck 
(P s ) durch Korrektur des reduzierten Querschnitts (A RRG ) durch eine KorrekturgroBe (AA RRG ) derart erfolgt, daB 
die Regelabweichung zwischen FuhrungsgroBe (P$_sen) und der ModellgroBe (P s ) minimiert wird. 

14. Verfahren nach Anspruch 12 und 13,, ' 
dadurch gekennzeichnet, 

daB die ModellgroBen vor dem Abgleich einer Mittelung (Gleichung 12. 1;1 3.1) unterzogen werden. 

15. Verfahren nach Anspruch 13, 
dadurch gekennzeichnet, 

daB dann, wenn die KorrekturgroBe (AA RRG ) einen vorgegebenen Schwellenwert in positiver Richtung uberschrei- 
tet, eine Erhohung der ModellgroBe Abgasgegendruck (P A ) erfolgt und wenn die KorrekturgroBe (AAj^ RG ) eine 
negative, applizierbare Schwelle unterschreitet, eine Erniedrigung der ModellgroBe Abgasgegendruck {P A ) erfolgt 

16. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daB die in den Zylinder einstrdmende Luftmasse (i,*,) durch die Beziehung 



berechnet wird, mit 



Abtastzeit oder Segmentzeit 

ModellgroBe des Luftmassenstromes wahrend des aktuellen Abtastschrittes oder Segments 
ModellgroBe des Luftmassenstromes wahrend des vergangenen Abtastschrittes oder Segmen- 

tes. 
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17. Verfahren nach Anspruch 1, * " 

dadurch gek nnz ichnet, 

daG die in d n Zylinder der Brehnkraftmaschine einstrdmende Frischgasmasse (mz^fAZ+H]) fur einen, bezuglich 
der aktuellen Lasterfassung zum Abtastzeitpunkt (N) in der Zukunft !iegenden : bestjmmten Vforhersagehorizont 
(H) auf der Grundlage der zeitlichen Anderungen der ModellgroBen Partialdrucke (P FG , P RG ) im Saugrohr (10) 
geschatzt wird. 


Claims 

1 . Method for determining the fresh air mass flowing into the cylinder of a combustion engine with 

an induction manifold (10) with a throttle valve (11) arranged within it, and a sensor device (14) which detects 
the degree of opening (DKW) of the throttle valve (19), 

an'exhaust gas recirculation device (1 9) with an exhaust gas recirculation valve (21 ) and a sensor device (22) 
which detects the degree of opening (AGRW) of the exhaust gas recirculation valve (21), 

a sensor (12, 13) which generates a load signal (rh DK LMM ; P S -sen) of tne combustion engine, 

an electronic control device which calculates a basic injection duration from the measured load signal 
(^dk.lmm^ p s_sen) and rotational speed (n) of the combustion engine, 

characterized in that 

the conditions in the induction manifold (10) are simulated by an induction manifold charging model (equations 
5.2 and 8,1), whereby at least the degree of opening (DKW) of the throttle valve (11), the degree of opening 
(AGRW) of the exhaust gas recirculation valve (21), the ambient pressure (Pj), the exhaust gas temperature 
(T A ), the temperature in the induction manifold (T s ) and the parameters representing the valve setting are 
used as input variables of the model, 

a model variable for the air mass flow (m DK ) at the throttle valve (11 ) and for the residual gas mass flow (^rc) 
at the exhaust gas recirculation valve (21) are formed in each case (equations 2.3 and 3.0) 
a model variable for the air mass flow into the cylinder (m z >i) or for the residual gas mass flow into the cylinder 
(^/.yi_Rc) are formed in each case as a function of the induction manifold pressure (P s ) or residual gas partial 
pressure (P RG ) (equations 4.0 and 7.0), 

The model variables for air mass flow at the throttle valve (rii DK ) : the residual gas mass flow at the exhaust 
gas recirculation valve (ih RC ) and air mass flow into the cylinder (rfu>i) are used as determining variables to 
calculate the induction manifold pressure (P s ) to determine the actual load of the combustion engine (equation 

5.2). 

the residual gas partial pressure (P RG ) in the induction manifold is calculated from model variables for the 
residual gas mass flow at the exhaust gas recirculation valve ( rh^) and residual gas mass flow into the 
cylinder (*z%ijtc) (equation 8,1 ), a ^ 
the fresh gas partjal pressure (P FG ) is determined from the induction manifold pressure (P s ) and residual gas 
partial pressure (P RG ) and 

the air mass (r^Fo) flowing A into the cylinder (17) is obtained by integration from the relationship between the 
fresh gas partial pressure (P FG ) and fresh air mass flow into the cylinder (rh* y ij.o). 

2. Method in accordance with claim 1 , 
characterized in that 

the model variables for the air mass flow (^dk) at the throttle valve (11) and for the residual gas mass flow (rn RC ) 
at the exhaust gas recirculation valve (21 ) are described by throttle settings with the aid of the flow-through equation 
for ideal gases. 

3. Method in accordance with claim 2, 
characterized in that 

the flow losses occurring at the throttling points (11,21) are allowed for in each case by using a reduced flow cross- 
section (A RDK , A RRG ) at the throttling points. 
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4. Method in accordance with claim 3, : - . • 

characterized in that 

the values for the reduced cross-sections (A RDK , A RRG ) are established on the engine test stand on the basis of 
steady-state measurements depending on the degree of opening (AGRW) of the exhaust gas recirculation valve, 
and stored in maps (KF1 , KF4) of a memory in the control device. 

5. Method in accordance with c\a\m 2, 
characterized in that 

in the representation of the model variables for the air mass flow (i^dk) at the throttle valve (11) and the model 
variable for the residual gas mass flow (^ RG ) at the exhaust gas recirculation valve (21), a flow-through function 
(Y) present in the flow-through equations is divided into single sections (i, j = 1 :..k) and these sections are ap- 
proximated by straight sections, whereby for the values of the gradient (mj, j), and for the absolute ^erm (n^ j) the 
particular straight sections, are determined as a junction of the ratio of induction manjfold pressure (P s ) to ambient 
pressure (Py) or induction manifold pressure (P s ) to exhaust gas back pressure (P/J and stored in a map in a 
memory in the control device. 

6. Method in accordance with claim 1, 
characterized in that 

the model variable jpr the air mass flow into the cylinder (*>*,) is described as a linear function of the induction 
manifold pressure (P s ) (equation 4.0) and gradient (y A ) and the absolute term (y 0 ) of the linear function is stipulated 
depending on at least one of the parameters, i.e. the rotational speed of the combustion engine, number of cylin- 
ders, induction manifold geometry, temperature of the air (T s ) in the induction manifold (10) or valve timing. 

7. Method in accordance with claim 1 , 
characterized in that 

the model variable for the residual gas mass flow in the cylinder (mz*_Rc) is described as a linear function of the 
residual gas partial pressure (P RG ) (equation 7.0) and the gradient (y n ) of the linear function is stipulated depending 
on at least one of the parameters, i.e. rotational speed of combustion engine, number of cylinders, induction man- 
ifold geometry, temperature of the air (T s ) in the induction manifold (10) or valve timing. 

B. Method in accordance with claim 6 or 7, 
characterized in that 

the parameters are determined by steady-state measurements on the engine test stand and stored in maps (KF2, 
• KF3) in a memory in the control device. 

9. Method in accordance with claim 1 , 
characterized in that 

the load signal (m DK _ LMM: Ps sen) measured by the load sensor (12, 13) is used in a model control circuit ^or 
correction and thus to adjust the air mass flow at the throttle valve (^dk) and induction manifold pressure (P s ) 
model variables; whereby the load signal (m DK LMM , P s _sen) serves as a reference variable for the control circuit. 

1 0. Method in accordance with claim 9, 
characterized in that 

the adjustment is carried out during steady-state and/or transient operation of the combustion engine and thus 
takes the response characteristic (equations 12.0 and 13.0) of the load sensor (12, 13) into account. 

11. Method in accordance with claim 10, 
characterized in that 

the response characteristic of the load sensor (12, 13) is modelled by a first order system with a working point- 
dependent delay time (T1 ). 

12. Method in accordance with claim 10, 
characterized in that 

a value of a reduced cross-section of the throttle valve (A RDK ) is assigned to each measured value of the throttle 
valve degree of opening (DKW) and the adjustment of the air mass flow at the throttle valve (rime) model variable 
is achieved by correction of the reduced cross-section (A RDK ) by a correction variable (AA RDK ) in such a way that 
the deviation between the reference variable (m DKJ _ MM ) and model variable (^ w ) is minimized. 
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13. Method in accordance with claim 10, • 
characterized in that 

a value of a reduced cross-section of the exhaust gas recirculation valve (A RRG ) is assigned to each measured 
value of the degree of opening (AGRW) of the exhaust gas recirculation valve (21) and the adjustment of the 
induction manifold pressure (P s ) model variable is achieved by correcting the reduced cross-section (A^ RG ) by 
means of a correction variable (AA RRG ) in such a way that the deviation between the reference variable (Ps-sen) 
and the model variable (P s ) is minimized. 

14. Method in accordance with claims 12 and 13, 
characterized in that 

the model variables are subjected to averaging (equations 12,1, 13.1) before the adjustment. . 

15. Method in accordance with claim 13, 

characterized in that if the correction variable (AA RRG ) exceeds a stipulated threshold value in the positive direction, 
an increase in the exhaust gas back pressure (P A ) model variable takes place and if the correction variable (AA RRG ) 
undershoots a negative, applicable threshold, a reduction in the exhaust gas back pressure (P A ) model variable 
results. 

1 6. Method in accordance with claim 1 , 
characterized in that 

the air mass (rh z „) flowing into the cylinder is calculated by the correlation 


m 2y ,[N] = (m Zyl [N - 1] + rfi 2yl [N]) 

where 

is the sampling interval or segment time 

is the model variable of the air mass flow during the actual sampling step or segment 

"* ^zyilN" 1 1 ia tne model variable of the air mass flow during the previous sampling step or segment 

17. Method in accordance with claim 1, 

characterized in that » 

the fresh gas mass (m Zy ,[N+H]) flowing into the cylinders of the combustion engine is estimated for a specific 
prediction horizon (H), relative to the actual load acquisition at sampling time point (N) in the future, on the basis 
of the change over time of the partial pressures (P FG , P RG ) model variables in the induction manifold (10). 

Revendications 

1. Procedede determination du debit massique d'airfrais introduit dans le cylindre d'un moteur a combustion interne, 
comprenant : 

un collecteur d'admission (10) et un papillon des gaz (11) dispose dans celui-ci, ainsi qu'un dispositif capteur 
(14) detectant le degre d'ouverture (DKW) du papillon des gaz (19), 

un dispositif de recyclage de gaz d'echappement (19) comportant une valve de recyclage de gaz d'echappe- 
ment (21), ainsi qu'un dispositif capteur (22) detectant le degre d'ouverture (AGRW) de la valve de recyclage' 
de gaz d'echappement (21), 

. - . un capteur (12; 13) produisant un signal de charge (m DK _ LMM ; P s _sen) du ^Xbux a combustion interne, . 
un dispositif de commande electronique qui calcule une duree d'injectidn de base & partir du signal de charge 
mesure (rn DK _ LMM ; P s _sen) et de ,a Vitesse de rotation (n) du moteur a combustion interne, 

caracterise en ce que : 


m Zy i[N] 
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- les conditions dans le collecteur d'admission (10) sont simu!6es au moyen cfun modele de tempi issage de 
collecteur d'admission (Equations 5.2, 8.1 ), au moins le degr§ d'ouverture (DKW) du papillon des gaz (11 ), le 
degre d'ouverture (AGRW) de la valve de recyclage de gaz d'echappement (21), la pression ambiante (P a ), 
la temperature de gaz d'echappement (T^), la temperature (T s ) dans le collecteur d'admission et des para- 
metres representant la position de valves etant utilises en tant que grandeurs d'entree du modele, " 
une grandeur de mod x I correspondant au d§bit massique d'air (m DK ) a I'endroit du papillon des gaz (11) et 
une grandeur de modele correspondant au debit massique de gaz residuels (a*) a I'endroit de la valve de 
recyclage de gaz d'echappement (21 ) sont respectivement formees (equations 2.3 ; 3.0), 
une grandeur de modele correspondant au debit massique d'air (a*n) dans le cylindre et une grandeur de 
modele correspondant au debit massique de gaz residuels (m 2yl m ) dans le cylindre sont respectivement for- 
mees en fonction de la pression de collecteur d'admission {P^ ou de la pression partielle de gaz residuels 
(P RG ) (equations 4.0 ; 7.0), 

a partir des grandeurs de modele que sont le debit massique (2« ) a I'endroit du papillon des gaz, le debit 
massique de gaz residuels (£„;) a I'endroit de la valve de recyclage de gaz d'echappement et le debit massique 
d'air (m 2yl ) dans le cylindre, sen/ant de grandeurs decisives pour la determination de la charge effective du 
moteur a combustion interne, la pression de collecteur. d'admission (P s ) est calculee (equation 5.2), 
a partir des grandeurs de modele que sont le debit massique de gaz residuels (a K )) a I'endroit de la valve de 
recyclage de gaz d'echappement et le debit massique de gaz residuels (Z Iyl K ) dans le cylindre, la pression 
partielle de gaz residuels {P RG ) dans ' e collecteur cTadmission est calculee (equation 8.1), A 
a partir de la pression de collecteur d'admission (P s ) et de la pression partielle de gaz residuels (P RG ). la 
pression partielle de gaz frais (Pp G ) est determined et 

a partir de la relation se presentant entre la pression partielle de gaz frais (P FG ) et le d6bit massique d'air frais 
(S, y :_«) dans le cylindre, le debit d'air (/tj^lfg) introduit dans le cylindre (17) est obtenu par integration. 

Procede suivant la revendication 1, caracterise en ce que les grandeurs de modele correspondant au d£bit mas- 
sique d'air (m 0K ) a I'endroit du papillon des gaz (11) et au debit massique de gaz residuels (m RC ) a I'endroit de la 
valve de recyclage de gaz d'echappement (21 ) sont decrites au moyen de Tequation d'ecoulement des gaz parfaits 
a travers des emplacements d'etranglement. 

Procede suivant la revendication 2, caracterise en ce qu'il est tenu compte de pertes de charge se presentant aux 
differents emplacements d'etranglement (11, 21) en utilisant des sections transversales d'ecoulement reduites 
(A RDK , A RRG ) respectives aux emplacements d'etranglement. 

Procede suivant la revendication 3, caracterise en ce que les valeurs correspondant aux sections transversales 
reduites (A RDK , A RFtG ) sont determinees sur banc d'essai de moteur a partir de mesures en regime permanent en 
fonction du degre d'ouverture (AGRW) de la valve de recyclage de gaz d'echappement et sont deposees dans 
des tables caracteristiques (KF1 ; KF4) d'une memoire du dispositif de commande. 

Procede suivant la revendication 2, caracterise en ce que, pour la representation des grandeurs de modele cor- 
respondant au debit massique d'air (m DK ) a I'endroit du papillon des gaz (11) et de la grandeur de modele corres- 
pondant au debit massique de gaz residuels a I'endroit de la valve de recyclage de gaz d'echappement (21), 
une fonction d'ecoulement (y) se presentant dans les equations d'ecoulement est divisee en differents segments 
(i, j = 1 . . . k) et ces segments font I'objet d'une approximation au moyen de segments de droite, les valeurs de la 
pente [n\ t j) et de I'abscisse a I'origine (nj j) des differents segments de droite^etant determinees sous forme de 
fonction respectivement du rapport de la pression de collecteur d'admission (P s ) a la pression ambiante (P u ) ou 
de la pression de collecteur d'admission (P s ) a la contre-pression de gaz d'echappement (P A ) et etant deposees 
dans une table caracteristique d'une memoire du dispositif de commande. 

Procede suivant la revendication 1 , caracterise en ce que la grandeur de modele correspondant au debit massique 
d'air (2 t/J ) dans le cylindre est decrite (equation 4.0) sous forme d'une fonction lineaire de la pression de collecteur 
d'admission (P s ), et la pente (y^) et I'abscisse a I'origine (y 0 ) de la fonction lineaire sont determinees en fonction 
d'au moins I'un des parametres que sont la Vitesse de rotation du moteur a combustion interne, le nombre de 
cylindres, les parametres geometriques de collecteur d'admission, la temperature (TS) de I'air dans le collecteur 
d'admission (10) et les durees de commande d'injecteur. 

Procede suivant ia revendication 1 , caracterise en ce que la grandeur de modele correspondant au debit massique 
de gaz residuels dans le cylindre est decrite (equation 7.0) sous forme d'une fonction lineaire de fa pression 

partielle de gaz residuels (P RG ), et la pente (7^ de la fonction lineaire sont determinees en fonction d'au moins 
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run des parametres que sont fa vilesse de rotation du moteur & combustion interne, le nombre de cylindres, les - 
parametres geometriques de collecteur d'admission, la temperature (TS) de Fair dans le collecteur d'admission 
(10) et les durees de commande d'injecteur. , 

s 8. Procede suivant la revendication 6 ou 7, caracterise en ce que les parametres sont determines au moyen de 
mesures en regime permanent sur banc d'essai de moteur et sont deposes dans des tables caracteristiques (KF2, 
KF3) d'une memoire du dispositif de commande. 

9. Procede suivant la revendication 1 , caracterise en ce que le signal de charge (^dkjlmm I P$_sen) mesurd par le 
10 capteur de charge (12 ; 13) sert, dans un circuit de regulation de modefe, k la correction et done k legalisation 

des grandeurs de modele que sont le debit massique d'air k I'endroit du papillon des gaz et la pression de 
collecteur d'admission (P s ), le signal de charge (m DK ^_ LMM ; P s _sen) servant de grandeur de pilotage du circuit de 
regulation. 

is 10. Procede suivant la revendication 9, caracterise en ce que legalisation est effectuee en regime permanent et/ou" 
en regime non permanent du moteur a combustion interne et il est alors tenu compte du comportement de trans- 
mission (equation 12.0 ; 13.0) du capteur de charge (12 ; ;13). 

1 1 . Procede suivant la revendication 1 0, caracterise en ce que le comportement de transmission du capteur de charge 
20 (12, 13) est modelise au moyen d'un systeme du premier ordre comportant un retard (T1) independant du point 

de travail. 

12. Procede suivant la revendication 10, caracterise en ce qu'une valeur respective d'une section transversale reduite 
(A RDK ) du papillon des gaz est associee a chaque valeur mesuree du degre d'ouverture de papillon des gaz (DKW) 

2S et 1'egalisation de la grandeur de modele qu'est le debit massique d'air im m )k I'endroit du papillon des gaz s'effectue 

en corrigeant la section transversale reduite (A m d au moyen d'une grandeur de correction (&A RDK ), de facon 
que I'ecart de regulation entre la grandeur de pilotage (mo/c.LMw) et la grandeur de modele (m Dh ) soit rendu minimal. 

1 3. Procede suivant la revendication 1 0, caracterise en ce qu'une valeur respective d'une section transversale reduite 
30 (A RRG ) de la valve de recyclage de gaz d'echappement est associee a chaque valeur mesuree du degre d'ouverture 

(AGRW) de la valve de recyclage de A gaz d'echappement (21) et I'egalisation de la grandeur de modele qu'est la 
pression de collecteur d'admission (P s ) s'effectue en corrigeant la section transversale reduite (A RRG ) au moyen 
d'une grandeur de correction (&A RRG ), de facon que I'ecart de regulation entre la grandeur de pilotage (P s _sen) 
et la grandeur de modele (P s ) soit rendu minimal. 

3S . 

14. Procede suivant la revendication 12 et 13, caracterise en ce que les grandeurs de modele sont soumises, avant 
I'egalisation, k une formation de moyenne (equation 12.1 ; 13.1). 

15. Procede suivant la revendication 13, caracterise en ce que, lorsque la grandeur de correction (AA RRG ) depasse 
40 vers le haut dans le sens positif une valeur de seuil prefixee, une augmentation de la grandeur de modele qu'est 

la contre-pression de gaz d'echappement (P A ) a lieu et, lorsque la grandeur de correction (&A RRG ) franchit vers 
le bas un seuil applicable negatif, une diminution de la grandeur de modele qu'est la contre-pression de gaz 
d'echappement (P A ) a lieu. 

45 16. Procede suivant la revendication 1 , caracterise en ce que le debit d'air (o^j.) introduit dans le cylindre est calcule 
au moyen de la relation 


so 


55 


m v [N] = & l)+m*[N]j 


dans laquelle 

TA est la duree d'analyse ou segment de temps, 

WzykW est ' a grandeur de modele du debit massique d'air pendant le pas d'analyse ou segment actuel, 

"Izy/W -1 ] est ,a grandeur de modele du debit massique d'air pendant le pas d'analyse ou segment passe. 
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1 7. Precede suivant la revendication 1 f caracterise en ce que le debit de gaz f rais (fffy/N+HI) introduit dans ie cylindre 
du moteur a combustion interne est estime, pour un horizon de prediction (H) determine situe dans le futur par 
rapport a la determination de charge actuelie a I'instant d'analyse^N), sur la base des variations dans le temps 
des grandeurs de modele que sont les pressions partielles (P FG , P RG ) dans le collecteur d'admission (10). 
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FIG 2 
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